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Abstract. We present an analysis of new and old high frequency VLBI data collected during the last 10 years at 5- 
22 GHz. For the jet components in the mas- VLBI jet, two component identifications are possible. One of them with 
quasi-stationary components oscillating about their mean positions. Another identification scheme, which formally 
gives the better expansion fit, yields motion with ~ 9 c for Ho = 65 kms" 1 Mpc -1 and go = 0.5. This model would 
be in better agreement with the observed rapid IDV and the expected high Lorentz-factor, deduced from IDV. 
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1. Introduction 

The S5 blazar 0716+714 is one of the most active BL Lac 
objects. It is extremely variable on time scales from hours 
to months. 0716+714 has yet no known redshift. Optical 
imaging of the underlying galaxy however reveals an esti- 
mate of z > 0.3 (Wagner et al. 1996). In the radio bands 
0716+71 is an intraday variable (IDV). It exhibits a very 
flat radio spectrum, extending up to at least 300 GHz. The 
variability appears to be correlated over wide ranges of the 
electromagnetic spectrum (Quirrenbach et al. 1991). The 
simultaneous variations between X-ray, optical and radio 
strongly suggest an intrinsic origin of variability. From 
IDV a high intrinsic Doppler factor could be expected. 
VLBI studies covering more than 20 years at 5 GHz show 
a core-dominated evolving jet extending to the north. The 
VLBI jet is oriented at 90° with respect to the VLA jet. 
At present it is not clear if 0716+714 is a superluminal 
source and how fast the motion of the VLBI components 
is. There are several possible models proposed with mo- 
tions between 0.05 — 0.5 mas yr -1 but some of them are 
based on relative unfrequently sampled data (Eckart et 
al. 1986, 1987; Witzel et al. 1988; Schalinski et al. 1992; 
Gabuzda et al. 1998; Tian et al. 2000). 

We present and discuss our preliminary results of 
the reanalysis of 10 years (25 epochs) of VLBI data on 
0716+714 at frequencies between 5 GHz and 22 GHz. Our 
analysis is based on four epochs from the CJF-Survey 
at 5 GHz between 1992 and 1999 QBritzcn et al. 1999Q , 
five epochs at 8.4 GHz from astrometric observations 



(|Ros et al. 2000|) and our own observations between 1994 
and 1999, five e pochs at 15 GHz from th e 2cm-Survey from 
1994 to 2001 flKcllcrmann et al. 199^ ) and 11 epochs at 
22 GHz from Jorstad et al. (2001) and our own data from 
1992 to 1997. 



2. Observations and results 

All the data were available as calibrated AIPS FITS files. 
The imaging of the source and the phase and amplitude 
self-calibration were done in Difmap. In Figure [|, we 
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Fig. 1. Top: Four epochs of 5 GHz maps. The maps are 
convolved with a 1.2 mas circular beam respectively. The 
lowest contours are of 2.5, 0.9, 2.8, 1.5 and 1.1 mJy /beam. 
Bottom: Five epochs 15 GHz maps. The maps are con- 
volved with a 0.5 mas circular beam respectively. The low- 
est contours are of 1.5, 2.5, 1.3 and 0.7mJy/beam. 



present a time sequence from the CJF-Survey and the 
2cm-Survey maps. The jet is clearly visible to the north at 
P.A. ~ 15°, but does not looks very straight. The ejection 
angle at the jet base shows regular variations and also the 
angle of the jet on larger scales seems to vary with time. 

2.1. Model-fits 

After imaging in Difmap we model fitted the self- 
calibrated data by circular gaussian components. We used 
circular components in order to reduce the degrees of free- 
dom and for simplify the analysis. In Figures |] & |], the 
core distance of individual VLBI components derived from 
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Fig. 2. Trajectories of model components in the quasi- 
stationary model. The lines are fits to illustrate the motion 
of identified components. 



our model-fits at various frequencies are plotted against 
time. A maximum error of 15% in the flux density can 
be estimated from the amplitude calibration and the un- 
certainties in the model-fits. The position error given in 
Figures ^| & [| is a superposition of the map cell size and 
an error of 10% of the core distance. 

3. Discussion 

3.1. Identification of the components 

To investigate the kinematics of 0716+714 we tried to 
identify individual model components along the jet by 
their properties, e.g. distance from the core, flux density 
and size. Since the observations were not phase-referenced, 
the absolute position information is lost, and it is impossi- 
ble to tell which component, if any, is stationary. Although 
we presently do not know if the brightest component is 
stationary we measured the positions with respect to the 
map phase center. But by a graphical analysis, which is 
presented in Figure || & |^, we could obtain two satisfactory 
scenarios for the kinematics of the jet. 

3.2. Scenarios 

Model 1: The first interpretation yields quasi-stationary 
components oscillating about their mean positions (see 
Figure |^) . This plot gives a reasonable argument presum- 
ing that the position of the brightest components is stable. 
If the core position were not stable all components at one 
epoch should have an offset with respect to the surround- 
ing observations which is not visible. 
Model 2: This model agrees with the theory that one can 
expect high Doppler factors for IDV sources. The com- 
ponents depart with O.lmasyr -1 , for the inner jet, and 
with up to 0.6 mas yr -1 , for the outer regions, from the 
core. Especially between 1993 and 1996, where we have 
many measurements spaced by only a few month, this 
model gives the most reasonable fit to the data. Assuming 




Fig. 3. Trajectories of model components in the superlu- 
minal motion scenario. The lines are fits to illustrate the 
motion of identified components. 

z = 0.3 these motions would correspond to a /3 app of 2-11 
(H = eSkms- 1 Mpc" 1 and q = 0.5). 

4. Conclusion 

At present we prefer Model 2 which would be in bet- 
ter agreement with the observed rapid IDV. Due to 
the observed radio optical correlation, the rapid IDV in 
0716+714 cannot be explained by refractive interstellar 
scattering. Therefore the source has to be very small 
(T B > 10 17 K) and the Doppler factor has to be > 50. The 
large oscillation of the components in Model 1 should make 
one suspicious about the identification. Further analysis 
with a comprehensive error calculation will show which 
model will prevail. 
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